Reentrant and Forward Phase Diagrams of 
the Anisotropic Three-Dimensional Ising Spin Glass 
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The spatially uniaxially anisotropic d = 3 Ising spin glass is solved exactly on a hierarchical lattice. 
Five different ordered phases, namely ferromagnetic, columnar, layered, antiferromagnetic, and spin- 
glass phases, are found in the global phase diagram. The spin-glass phase is more extensive when 
randomness is introduced within the planes than when it is introduced in lines along one direction. 
Phase diagram cross-sections, with no Nishimori symmetry, with Nishimori symmetry lines, or 
entirely imbedded into Nishimori symmetry, are studied. The boundary between the ferromagnetic 
and spin-glass phases can be either reentrant or forward, that is either receding from or penetrating 
into the spin-glass phase, as temperature is lowered. However, this boundary is always reentrant 
when the multicritical point terminating it is on the Nishimori symmetry line. 
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I. INTRODUCTION 

The Ising spin glass [l[ yields a phase diagram with 
a distinctively complex ordered phase, in d = 3. A 
wide accumulation of methods and results has occurred 
for this system. Most remarkably, in spite of its high 
spatial dimension and complex ordering behavior, ex- 
act or precise information is being obtained for this 
system, p, S, 1, |, 0, 0, S, S, M, EL d Thus, in the 
phase diagram in terms of temperature and concentra- 
tion of antiferromagnetic bonds, the occurrence of the 
Nishimori symmetry line has been deduced @, Q and the 
accurate location of the multicritical point has been pre- 
dicted [HJII2]- Furthermore, in systems with the Nishi- 
mori symmetry, it has been shown that the ferromagnetic 
phase cannot extend to antiferromagnetic bond concen- 
trations beyond that of the multicritical point. @, Q The 
two remaining options being a strai ght line or a reen- 
trance situation, subsequent works [l9t |20| on hierar- 
chical lattices have shown that for these systems, the 
spin-glass phase diagram is reentrant, namely that be- 
low the multicritical point, the ferromagnetic phase re- 
cedes from the spin-glass phase as temperature is low- 
ered. Exact results recently have also been extended to 
Potts spin glasses. [2lj These results complement recent 
precise calculations, using Monte Carlo simulations, on 
cubic lattices. 0, Q M, M, E3, 01 

A spatially uniaxially anisotropic d = 3 system is 
studied in this work, to our knowledge the first study 
of quenched randomness and frustration in a spatially 
anisotropic higher-dimensional system. In fact, both 
anisotropy and quenched randomness have acquired in- 
creased relevance from high-temperature superconduc- 
tivity results. [22], [H| Our calculation is exact for a hier- 
archical lattice and approximate for a cubic lattice. We 
find a rich phase diagram (e.g., FigJT]) with five differ- 
ent ordered phases, namely with ferromagnetic, antifer- 



romagnetic, layered, columnar, and spin-glass order. The 
spin-glass phase is more extensive when randomness is in- 
troduced within the planes than when it is introduced in 
lines along one direction. 

The global phase diagram includes cross-sections with 
no Nishimori symmetry, cross-sections with Nishimori 
symmetry lines, and a cross-section entirely imbedded 
within Nishimori symmetry. Thus, the multicritical point 
between the spin-glass, ferromagnetic, and paramagnetic 
phases, previously found to occur on the Nishimori sym- 
metry line, is also found here at points with no Nishi- 
mori symmetry, but renormalizes to a fixed distribution 
of interaction probabilities that obeys Nishimori symme- 
try. Nevertheless, we find that the boundary between the 
ferromagnetic and spin-glass phases can be either reen- 
trant or forward, that is either receding from or penetrat- 
ing into the spin-glass phase, as temperature is lowered. 
When the multicritical point is not on the Nishimori sym- 
metry line, the ferromagnetic-spin glass boundary can be 
reentrant or forward. However, when the multicritical 
point is on the Nishimori symmetry line, this boundary 
is always reentrant [19l.l20l|. consistently with the rigorous 
result 



II. UNIAXIALLY ANISOTROPIC SPIN GLASS 

The uniaxially anisotropic Ising spin-glass system has 
the Hamiltonian 



u (ijU 



(1) 



where Sj = ±1 at each site i, (ij) u denotes a sum over 
nearest-neighbor pairs of sites along the z direction (u — 
z) or in the xy plane (u = xy), and the bond strengths 
Kfj are equal to K u > with probability 1 — p u and 
—K u with probability p u , respectively corresponding to 
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FIG. 1: (Color on-line) Constant-temperature cross-sections of the global phase diagram for K z /K xy = 0.5, as a function 
of p xy and p z , which are the concentrations of antiferromagnetic xy and z bonds, respectively. At low temperatures (high 
K xy ), the central spin-glass (SG) phase separates the corner ferromagnetic (F), columnar (C), antiferromagnetic (A), and 
layered (L) phases. The diagrams are twofold symmetric along each axis, but not fourfold symmetric, due to the difference 
between longitudinal (p xy = 0) and transverse (p z = 0) spin glasses. As temperature increases, the paramagnetic (P) phase 
appears at the central point, first reaches the transverse spin-glass system and eliminates the spin-glass phase, then reaches 
the longitudinal spin-glass system and eliminates the spin-glass phase. In the latter system, the spin-glass and paramagnetic 
phases simultaneously occur for a very narrow range of temperatures, as also seen in the inset in the lower left panel of Fig. 3. 
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FIG. 2: Construction of the uniaxially anisotropic d=3 hier- 
archical model. Two graphs are mutually and repeatedly self- 
imbedded. Note that for K xy = 0, K z = 0, and K xy = K z , 
the system reduces respectively to the d = 1, isotropic d — 2 
and d = 3 systems. 



riety of problems, including chaotic rescaling [23, l28j . 
spin-glass [l9j . random- field [29| . Schrodinger equation 
[3fj| , lattice- vibration [3l[ , dynamic scaling |32| , random- 
resistor network (33|, aperiodic ma gne t [34| , complex 
phase diagram [35[ , directed-path [361.1371] , heteropolymer 
HH, directed-polymer [39j | , an d, most recently, scale-free 
and small-world network gO, El, E HI H ES EE E3 
systems, etc. More recently, hierarchical lattices have 
been created [48[ for the study of spatially anisotropic 
systems. The mutual repeated self-imbedding of two 
appropriately chosen graphs, with differentiated interac- 
tions, yields a uniaxially anisotropic system, whereas a 
higher number of g raphs is needed to achieve higher spa- 
tial anisotropy. *48j These hierarchical systems must re- 
duce to isotropy and/or lower spatial dimensions when 
corresponding interactions are set equal to each other or 
to zero, as illustrated in Fig[2] An anisotropic hierar- 
chical lattice has already been used to obtain the phase 
diagram of the uniaxially anisotropic d = 3 tJ model of 
electronic conduction. [22| When imbedded into the hi- 
erarchical lattice of FigJ^l the Hamiltonian ([T]) yields a 
uniaxially anisotropic d — 3 spin-glass system that is ex- 
actly soluble. 



ferromagnetic and antiferromagnetic interaction. When 
imbedded into a cubic lattice, the Hamiltonian ([I]) yields 
a uniaxially anisotropic d — 3 system. 

Hierarchical lattices are d-dimensional lattices yielding 
exact renormalization-group solutions to complex statis- 
tical mechanics problems. These lattices are constructed 
by the repeated self-imbedding of a graph into a bond 
[HlHi, 26]. The shortest path between the external ver- 
tices of the graph gives the length rescaling factor b and 
the number of bonds in the graph gives the volume rescal- 
ing factor b d , from which the dimension d is determined. 
Hierarchical lattices have been used to study a wide va- 



III. EXACT RENORMALIZATION-GROUP 
SOLUTION: FLOWS OF THE QUENCHED 
DISTRIBUTIONS OF THE ANISOTROPIC 
SPIN-GLASS INTERACTIONS 

The renormalization-group solution proceeds in the 
direction opposite to the construction of a hierarchical 
model. Each graph is replaced by a renormalized bond 
via summation over the spins on the internal sites of the 
graph. This is achieved by a combination of two types 
of steps: the replacement, by a single bond K^j, of two 
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bonds that are either in parallel, referred to as bond- 
moving: 

K ij =K! i + K l i f, (2) 
or in series, referred to as decimation: 

~ i r cosh(i^ + if Jfc ) i 

ife 2 [ca&HKij - K ik )\ ' y ' 

The quenched probability distribution V(K) of the re- 
placing bond is calculated by the convolution 

V{K) = J dK I dK II Vi(K I )Vii(K II )S(K-R(K I ,K 11 )), 

(4) 

where R(K ] K H ) is the right-hand side of Eq.@ or ©, 
K 1 and K 11 are the interactions entering the right-hand 
side of either of these equations, with quenched proba- 
bility distributions Vi{K r ) and V U {K U )^, [n| 

Accordingly, the renormalization of V xy is obtained as 
follows, following the upper FigTJJin the direction oppo- 
site to the arrow: (i) from the bond-moving of V xy with 

itself, obtaining V\ ; (ii) from the bond- moving of V z with 
itself, obtaining V% (hi) from the decimation of V\ and 
■Pi, obtaining V3] (iv) from the decimation of Vi and Vi, 
obtaining Vn\ (v) from the decimation of V xy and V xy , 
obtaining V5] (vi) from the decimation of V3 and V\ 1 
obtaining Ve ; (vii) from the decimation of Va and V2, 
obtaining V7; (viii) from the decimation of V5 and V xy , 
obtaining Vs', (ix) from the bond-moving of V§ and V7, 
obtaining Vg; (x) from the bond-moving of V7 and V7, 
obtaining V10; (xi) from the bond-moving of Vg and V10, 
obtaining V\\ (xii) finally, from the bond-moving of V\\ 
and Vs, obtaining the renormalized quenched distribu- 
tion V' xy . Thus, in each renormalization-group step, the 
renormalized distribution V' xy is obtained from the con- 
volutions of 27 unrenormalized distributions V xy and V z . 
The renormalized distribution V' z is similarly obtained 
from the convolutions of 27 unrenormalized distributions 
V xy and V z , but with a different sequencing dictated by 
the lower FigJ^l 

The renormalization-group transformations of the 
quenched probability distributions V xy and V z , given 
in the preceding paragraph, are implemented numeri- 
cally, resulting in a distribution of interaction-strength 
values and a probability associated with each value, 
namely a histogram. Thus, the initial ±K U double- 
delta distribution functions, described after Eq.(fT]), are 
of course not conserved under the scale coarsening of 
the renormalization-group transformation. The number 
of histograms increases after each convolution. When a 
maximum number of histogra ms, set by us, is reached, a 
binning procedure is applied [l^, H^] : Before each con- 
volution, the range of interaction values is divided into 
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TABLE I: Sinks of the renormalization-group flows in the dif- 
ferent phases. These sinks are characterized here in terms of 
the average positive and negative interactions of their limiting 
quenched probability distribution. 



bins, separately for positive and negative interactions. 
The interactions falling into the same bin are combined 
according to their relative probabilities. The convolution 
then restores the set maximum number of histograms. In 
this work, we have used the maximum number of 90,000 
for histograms for each distribution V xy and V z . 



IV. PHASE DIAGRAMS AND FIXED 
DISTRIBUTIONS 

We have obtained the global phase diagram of the uni- 
axially anisotropic d = 3 spin-glass system in terms of the 
original interactions and probabilities (K xy , K z ,p xy ,p z ). 
In each thermodynamic phase, quenched probability dis- 
tributions flow, under repeated renormalization-group 
transformations, to a limiting behavior (sink) charac- 
teristic of that thermodynamic phase. Phase bound- 
ary points flow to their own characteristic (unstable) 
fixed distributions, shown below. Analysis at these un- 
stable fixed distributions yields the order of the phase 
transitions. [29| 

We find six different phases for this system, with cor- 
responding sinks characterized in Table I in terms of the 
average positive and negative interactions of the limiting 
distribution. These phases are the ferromagnetic, antifer- 
romagnetic, layered, columnar, spin-glass ordered phases 
and the disordered paramagnetic phase. In the layered 
phase, the spins are mutually aligned in each xy plane; 
these planes of mutually aligned spins form an antifer- 
romagnetic pattern along the z direction. In the colum- 
nar phase, the spins are mutually aligned along the z 
direction; these lines of mutually aligned spins form an 
antiferromagnetic pattern along the xy directions. Both 
of these phases are thus distinct from the antiferromag- 
netic phase, which is antiferromagnetic in all three di- 
rections. There is a single spin-glass phase, extending to 
anisotropic systems. 



A. Phase Diagrams with no Nishimori Symmetry 

Cross-sections of the global phase diagram are given 
in Figs[U 121 and 2J All phase transitions in these fig- 
ures are second order. FigfTJ shows constant-temperature 
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FIG. 3: (Color on-line) Temperature-concentration phase di- 
agrams for isotropically mixed (upper left), transverse (upper 
right), longitudinal (lower left), and p xy = 0.5p z spin-glass 
systems. In all cases, K z /K xy — 0.5. The upper left and 
right phase diagrams are seen to be, respectively, reentrant 
and forward, namely with a ferromagnetic phase that, respec- 
tively, recedes from or proceeds towards the spin-glass phase 
as temperature is lowered, as clearly seen in the insets. There 
are no points obeying Nishimori symmetry in the phase dia- 
grams of this figure. Note the remarkably narrow spin-glass 
phase, reaching zero-temperature, in the longitudinal spin- 
glass system, as also seen in the inset. All phase transitions 
in this figure are second order. 



1 z 1 xy 
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FIG. 4: (Color on-line) Zero-temperature phase diagrams of 
the longitudinal (left column) and transverse (right column) 
spin-glass systems. With the appropriate reversal in variables, 
the transverse and longitudinal spin-glass phase diagrams are 
seen here to be qualitatively similar, but quantitatively differ- 
ent. The spin-glass phase is more extensive in the transverse 
case. All phase transitions in this figure are second order. 



cross-sections of the global phase diagram as a func- 
tion of p xy and p z . At low temperatures (high K xy ), 
the central spin-glass (SG) phase separates the corner 
ferromagnetic (F), columnar (C), antiferromagnetic (A), 
and layered (L) phases. The diagrams are twofold sym- 
metric along each axis, but not fourfold symmetric, due 
to the difference between transverse (p z — 0) and lon- 
gitudinal (p X y = 0) spin glasses. As temperature in- 
creases, the paramagnetic (P) phase appears at the cen- 
tral point, first reaches the transverse spin-glass system 
and eliminates the spin-glass phase, then reaches the lon- 
gitudinal spin-glass system and eliminates the spin-glass 
phase. Fig[3] shows temperature-concentration phase di- 
agrams for isotropically mixed, transverse, longitudinal, 
and p xy — 0.5p z spin-glass systems. The upper left and 
right phase diagrams are seen to be, respectively, reen- 
trant and forward, namely with a ferromagnetic phase 
that, respectively, recedes from or proceeds towards the 
spin-glass phase as temperature is lowered, as clearly 
seen in the insets. The Nishimori symmetry (see below) 
is obeyed only at four isolated ordinary points in each 
cross-section in Fig[T] and is not obeyed at any point in 
the phase diagrams in Figs. [3J [H so that the forward 
behavior is not excluded by the rigorous results 0, Q . 

A remarkably narrow spin-glass phase, reaching zero- 
temperature, occurs in the longitudinal spin-glass sys- 
tem. Zero-temperature phase diagrams are shown in 
Fig|4] for the longitudinal (left column) and transverse 
(right column) spin-glass systems. With the appropriate 
reversal in variables, the longitudinal and transverse spin- 
glass phase diagrams are seen in this figure to be qual- 
itatively similar, but quantitatively different. The spin- 
glass phase is more extensive in the transverse case. This 
can be understood from the more extensive intermixing 
of the ferromagnetic and antiferromagnetic bonds. 



B. Temperature-Concentration Phase Diagrams 
with Nishimori Symmetry Curved Lines 

The Nishimori symmetry condition 0, 0] for isotropic 
systems, 



1-P 

P 



(5) 



generalizes, for uniaxially anisotropic spin-glass systems, 
to 



L_P*!t = e ±2*., and L^£ 



( 6 ) 

Pxy Pz 

For Nishimori symmetry to obtain, both equations have 
to be satisfied, but the signs in the exponents can be 
chosen independently. The Nishimori condition, in its 
general form 



Vu{-K u ) = ±2Ku 
Vu(Ku) 



(7) 
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FIG. 5: (Color on-line) Phase diagrams with Nishimori sym- 
metry lines (dashed) for different anisotropy parameters: The 
ratio K z /K xy is 2, 1 and 0.5 from top to bottom. In the left 
column, p xy satisfies the Nishimori condition. In the right 
column, K z satisfies the Nishimori condition. All phase tran- 
sitions in this figure are second order. 

for each histogram pair of each distribution, is invari- 
ant (closed) under our renormalization-group transfor- 
mation. 

If one of the two conditions in Eq.((6|) is fixed, phase 
diagram cross-sections are obtained, in which Nishimori 
symmetry holds along a line. Thus, throughout the three 
phase diagrams on the left in Fig. 5, the condition on 
(K xy ,p xy ) is fixed. The condition on (K z ,p z ), and there- 
fore Nishimori symmetry, is satisfied along the dashed 
lines on the left in Fig. 5. In these temperature versus 
concentration phase diagrams, it is seen that the multi- 
critical points between the ferromagnetic, spin-glass, and 
paramagnetic phases lie on the Nishimori symmetry line. 
Furthermore, it has been proven 0, [1[ that a forward 
phase diagram cannot occur below such a multicritical 
point that is on the symmetry line. On the left in Fig. 5, 
this is indeed the case, with reentrant phase diagrams, as 
also seen in isotropic spin glasses [III H3| • Recall that in 
Sec.IVA, multicritical points, between the same phases as 
here, that do not lie on Nishimori symmetry occur with 
both reentrant and forward phase diagrams. However, 
the latter non-symmetric multicritical points flow, un- 
der renormalization-group transformations, to the (dou- 
bly unstable) fixed distribution of the symmetric multi- 
critical points, therefore being in the same universality 
class and having the same critical exponents. 

In the three phase diagrams on the right of Fig. 5, the 
condition on (K z ,p z ) is fixed. In these concentration- 
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FIG. 6: (Color on-line) The Nishimori condition for K z is 
held throughout the leftmost figure and for K xy throughout 
the center figure. The complimentary Nishimori condition, 
for K xy and K z respectively, is held along the dashed straight 
lines, which intersect the ordered (F, L, A, or C)-spinglass- 
paramagnetic multicritical points. In the rightmost figure 
both conditions are satisfied throughout the figure. In this 
figure, the phase boundaries around the paramagnetic phase 
are actually lines of the multicritical points where the para- 
magnetic, ordered (F, L, A, or C), and spin-glass (not seen 
in this cross-section) phases meet. In the side figures, first- 
order boundaries (dotted) occur between the ferromagnetic 
and layered phases, and between the antiferromagnetic and 
columnar phases, terminating at d — 2 critical points. All 
other phase transitions (full lines) in this figure are second 
order. 



concentration phase diagrams, the multicritical points 
between the ordered (ferromagnetic, antiferromagnetic, 
layered, or columnar), spin-glass, and paramagnetic 
phases again lie on the Nishimori symmetry lines. 



C. Concentration-Concentration Phase Diagrams 
with Nishimori Symmetry Straight Lines 

In the phase diagrams in Fig. 5, the ratio K z /K xy is 
held constant. On the left and center of Fig. 6, again the 
condition in Eq.© on one interaction is fixed and the 
other interaction strength is held constant. Thus, the 
Nishimori symmetry lines becomes straight lines. The 
multicritical points between the ordered (ferromagnetic, 
antiferromagnetic, layered, or columnar), spin-glass, and 
paramagnetic phases again lie on the Nishimori symme- 
try lines. In the left phase diagram, due to the en- 
forced Nishimori symmetry condition, K z = along 
the line p z = 0.5 and the system reduces to d = 2. 
Along this line, first-order transitions between ferromag- 
netic and layered phases and between antiferromagnetic 
and columnar phases terminate at d = 2 critical points. 
From p z ^ 0.5, d — 3 second-order boundaries between 
each ordered phase and the paramagnetic phase termi- 
nate on the d — 2 critical points. In the center phase 
diagram, due to the enforced Nishimori symmetry con- 
dition, K xy — along the line p xy = 0.5 and the system 
reduces to d = 1. Accordingly, the system is disordered 
(paramagnetic) along the entire length of this line. 
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FIG. 7: (Color on-line) Fixed distributions, with circles 
and crosses showing one renormalization-group transforma- 
tion and thereby by their exact superposition attesting to the 
fixed nature of the distributions. The distributions have been 
binned for exhibition purposes, (a) For the ferromagnetic- 
spinglass phase boundary, a runaway to infinite coupling; 
(b) for the paramagnetic-spinglass phase boundary. Both 
of these fixed distributions are spatially isotropic, attract- 
ing isotropic and anisotropic boundaries, and do not obey 
Nishimori symmetry, (c) For the ferromagnetic-spinglass- 
paramagnetic multicritical point. This fixed distribution 
is spatially isotropic and obeys Nishimori symmetry. This 
fixed distribution attracts the isotropic multicritical point, 
which obeys Nishimori symmetry, and anisotropic multicriti- 
cal points, which obey and do not obey Nishimori symmetry. 
The fixed distributions for the antiferromagnetic-spinglass, 
columnar-spinglass, layered- spinglass phase boundaries and 
for the antiferromagnetic-spinglass-paramagnetic, columnar- 
spingiass-paramagnetic, layered-spinglass-paramagnetic mul- 
ticritical points are as shown here in (a) and (c) respectively, 
but with the appropriate K xy — > —K xy and/or K z — > — K z 
reflections, (d) Fixed distribution for the spin-glass phase. 
This phase sink is an isotropic runaway, attracting both spa- 
tially isotropic and anisotropic spin-glass phase points, and 
does not obey Nishimori symmetry. 



D. The Phase Diagram Entirely Imbedded in 
Nishimori Symmetry 

In the rightmost Fig. 6, both conditions of Eq.© are 
satisfied throughout the figure. With two symmetry con- 
straints, this is a unique surface in the global phase di- 
agram of our model. The system reduces to d — 2 and 
d = 1, as explained above, for p z — 0.5 and p xy = 0.5 
respectively. The phase boundaries around the paramag- 
netic phases are actually lines of the multicritical points 
where the paramagnetic, ordered (ferromagnetic, layered, 
antifcrromagnetic, or columnar), and spin-glass (not seen 



in this cross-section) phases meet. 

No spin-glass phase occurs within the Nishimori- 
symmetric subspace. The phase transitions seen 
in the rightmost Fig. 6, namely ordered-spinglass- 
paramagnetic multicritical and ferromagnetic-layered, 
antiferromagnetic-layered first-order transitions, are the 
only phase transitions of the system that occur under 
Nishimori symmetry. 



E. Fixed Distributions 

The fixed distributions underpinning the phase 
diagrams of this system are given in Fig. 7. The 
fixed distributions for the ferromagnetic-spinglass 
boundary, paramagnetic-spinglass boundary, and 
the ferromagnctic-spinglass-paramagnetic multicrit- 
ical points are spatially isotropic, but attract both 
spatially isotropic and anisotropic phase transitions. 
The fixed distribution for the ferromagnetic-spinglass- 
paramagnetic multicritical points obeys Nishimori 
symmetry, but attracts multicritical points that obey 
and do not obey Nishimori symmetry In the latter 
cases, as seen above, both reentrant and forward phase 
diagrams occur. 

The fixed distributions for the antiferromagnetic- 
spinglass, columnar-spinglass, layered-spinglass phase 
boundaries and for the antiferromagnetic-spinglass- 
paramagnetic, columnar-spinglass-paramagnetic, 
layered-spinglass-paramagnetic multicritical points 
are as shown in Fig. 7 (a) and (c) respectively, but 
with the appropriate K xy — > —K xy and/or K z — > — K z 
reflections. 



V. CONCLUSION 

The exact solution of the spatially uniaxially 
anisotropic spin glass on a d = 3 hierarchical lattice 
yields new phase diagrams. In view of the semiquan- 
titative agreement between spatially isotropic spin-glass 
results on cubic and hierarchical lattices [19(, it would 
certainly be worthwhile to investigate on cubic lattices 
the new phenomena found in the present study. Further- 
more, the exact study of spin glasses on fully anisotropic 
d — 3 hierarchical lattices [48| may yield even more new 
phase transition phenomena. 
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